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ABSTRACT 

The discovery of (3 Cephei stars in low metallicity environments, as well as the dif- 
ficulty in theoretically explaining the excitation of the pulsation modes observed in 
some (3 Cephei and hybrid SPB-/3 Cephei pulsators, suggest that the "iron opacity 
bump" provided by stellar models could be underestimated. We analyze the effect 
of uncertainties in the opacity computations and in the solar metal mixture, on the 
excitation of pulsation modes in B-type stars. We carry out a pulsational stability 
analysis for four grids of main-sequence models with masses between 2.5 and 12 Mq 
computed with OPAL and OP opacity tables and two different metal mixtures. 

We find that in a typical (3 Cephei model the OP opacity is 25% larger than OPAL 
in the region where the driving of pulsation modes occurs. Furthermore, the difference 
in the Fe mass fraction between the two metal mixtures considered is of the order 
of 20%. The implication on the excitation of pulsation modes is non-negligible: the 
blue border of the SPB instability strip is displaced at higher effective temperatures, 
leading to a larger number of models being hybrid SPB-/3 Cephei pulsators. Moreover, 
higher overtone p-modes are excited in [3 Cephei models and unstable modes are found 
in a larger number of models for lower metallicities, in particular (3 Cephei pulsations 
are also found in models with Z=0.01. 

Key words: radiative transfer - stars:abundances - stars:early-type - stars:interiors 
- stars:oscillations - stars:variables:other 



1 INTRODUCTION 

At the beginning of th e nineties, OPAL opacity tables 
iRoeers fe Igfesiaa Il992l) signified a revolution in stellar 
physics. In particular, an increase of the Rosseland mean 
opacity (kr) as large as a factor 3 for temperatures 
near 3 10 5 K (known as "Z-bump") allowed to solve 
the long-standing problem of B-type pulsators: /3 Cep 
and Slowly Pulsating B-stars (SPB) pulsate due to k- 
mechanism activated by this metal opacity bump l|Cox_et_alJ 
19921: iKiriakidis et alJll992t iMoskalik fe Dziembowskilll99a 
Dziembowski et al.l 1199311 . Subsequent improvements in 
opacity calculations by including other iron-group elements 
in the metal mixture and intermediate-coupling for Fe tran- 
sitions, have led to an additional enhancement of opacity 
at logT ~ 5.2 and low densities, and therefore, to a de- 
crease of the still remaining discrepancies between theoreti- 
cal pulsation models and observa tions. The stability analysis 
of B-stars bv lPamvatnvkhl lll99 9) for models computed with 
updated OP AL and OP opacity tables (Opacity Project, 
ISeatonlll996t and references therein) allowed to conclude 



that theoretical standard models were able to explain most 
of the observed ft Cep and SPB stars. 

However, as observational capabilities progress, B-type 
pulsat ors are now found in low metallicity environments (see 
e.g. iKolaczkowski et alJl20od . and references therein) , and a 
large number of pulsation modes is being detected in B stars, 
with frequency domains sometimes revealing new discrepan- 
cies between theory and observations. For instance, the two 
(3 Cep stars 12 Lacertae and v Eridani present low order p- 
modes with frequencies larger than those predicted by pulsa- 
tion models, as well as high-order g-modes (SPB typ e oscilla- 
tion) iHandler et al 1l200rj:ljerzvk iewicz et al . 2005, and ref- 
erences therein). In order to expla in these pulsation features 
in v Eri. lPamvatnvkh et al] ]2004]l propose d a local enhance - 
ment of iron in the "Z-bump" region, like lCox et al] <ll992l) 
suggested to overcome the difficulties to explain f3 Cep pul- 
sations with the first OPAL opacity tables. 

Regardless of whether other additional physical pro- 
cesses must be included in the B-type stellar modeling, we 
propose here to first analyze the uncertainties in the predic- 
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Figure 1. Comparison between opacities computed with OPAL, OP and with different metal mixtures in the structure of a 4 and 10 
Mq models on the ZAMS. 



tions of B-type pulsation models due to the uncertainties in 
the basic physical inputs i n the standard model. 

It is known (see e.g. iMoskalik fc DziembowskH I1992T) 
that stellar pulsation models are sensitive to ~20% changes 
in opacity. Recent improvements in OP opacity computa- 
tions have, on the one hand, increased the opacity in the 
regions where OP values were previously much lower than 
OPAL values, leading to a good agreement between both 
calculations. On the other hand, the new atomic data pro- 
vide a 18% enhancement of opacity in the Z-bump region. 

Besides, the chemical composition of the Sun and of B- 
ty pe stars is still matte r of lively debate (see, e.g. the works 
of Usplund et al.ll2005l and lCunha et alJl2006lL In this work 
we therefore investigate the impact that the uncertainties 
about the metal mixture have on the instability strips of 
B-type stars, as well as on the frequency domain of excited 
modes. 



2 METAL MIXTURES 

The r ecent re-analysis of solar spectrum bv lAsplund et ail 
J2005I) (thereafter AGS05) including NLTE effect as well as 
tri-dimensional model atmosphere computations, has led to 
a significant decrease of C, N, O and Ne solar abundances 
leading to a solar metallicity 30% smaller than the value pro- 
vided by the "standard" GN93 mixture IIGrevesse fc Noelsl 
Il993t) . The corresponding decrease in solar opacity (~ 20%) 
at the base of the convective envelope ruins the good agree- 
ment between the standard solar model and the seismic 
one. These new CNO abundances, however, agree with spec- 
troscopic measurements in B-type stars, and solve the old 
discrepancy between a supermetallic Sun and the chemical 
composition of B-type stars and HII r egions in the solar 
environment (iTurck-Chieze et al-lliool ICunha et alJl200rl 
and references therein). On the other hand, while the lower 



O abundance in AGS05 mixture decreases the opacity at 
the bottom of the convective envelope and increases the dis- 
crepancy between standard solar model and helioseismology, 
an iron mass fraction 25% larger in the new mixture with 
respect to GN93 one, will favorably affect the excitation of 
f3 Cep and SPB pulsation modes in early-type stars. 

Among the different solutions suggested to recover the 
agreement between standard model and helioseismology, an 
enhancemen t of the Ne abundance by a factor ~3.5 with re- 
sp ect to thelAsplund et alJ 112005!) value ha s been proposed 
bv lAntia fc Basul J2005h ~and iBahcall et alJ fcOOElL The in- 
crease of Ne would compensate for the opacity lost due 
to the O abundance drop. The problem is th at Ne cannot 
be dir ectly measured in the solar photosphere. ICunha et alJ 
(2006) suggested that determining non-LTE Ne abundances 
in a sample of eleven B-stars members of the Orion asso- 
ciation could help solving the problem of Ne abundance. 
The studied stars span an effective temperature range from 
20000 K to 29000 K. The result of this wor k is a Ne abun- 
dance in B-type stars 0.3 dex larger than in lAsplund et all 
(2005) mixture, but CNO in agreement with the new solar 
abundances. 

In the following sections we will adopt [ Asp lund et alJ 

d2005f) mixture with A(Ne)=8.11 from lCunha et alJ (1200^) 
as the new metal mixture (AGS05+Ne) for B-type stars, 
and we will compare results obtained with the "standard" 
GN93. The re-normalization of AGS05+Ne to keep the same 
molecular weight results finally in an enhancement of ~ 20% 
in the iron and nickel mass fractions with respect to GN93. 
In Fig. the dotted curves show the relative differences in 
the Rosseland mean opacity between GN93 and AGS05+Ne 
(for the same metal mass fraction Z=0.02) for the structure 
of a 4 and a 10 Mq models. The effect of the enhancement 
of Fe and Ni is to increase the height and the width of the 
Z-bump of the Rosseland mean opacity (kr). The conse- 
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quences on the excitation of pulsation models are described 
in Sec.|S] 



3 OPACITIES: OPAL VERSUS OP 

OP opacity tables have been recently updated including 
data for the inner-shell transitions, inter-combinati on lines 
and improvements in pho toionization cross-sections llSeatonl 
120051: iBadnell et all 12005). These changes result in an en- 
hancement of kr at high density and temperatures, and an 
increase of 18% of opacity in the Z-bump due to the new 
Fe atomic data. The new OP opacity tables ar e much closer 
to th e OPAL ones than they were previously llSeaton et alJ 
Il994l) . In particular, the differences at high temperature and 
high density have almost disappeared. Some differences re- 
main, however, at low temperatures (logT < 5.5): the OP 
Z-bump in kr presents a hot wing slightly larger than the 
OPA L one. The com parisons between OP and OPAL by 
IBadnell et all J2005I) at fixed log R (R = p/T 6 3 ) show that 
for high densities the differences are of the order of 5-10%, 
but at low density (logi? ~ ~4), differences up to 30% ap- 
pear due to the differences in the high-temperature wing of 

the Z-bump. 

] ' r previous OP tables dSeaton et alJ Il994h already 
presented a Z-bump in kr shifted by 15000-20000 K to 
higher temperatures compared with OPAL. The effect of 
this difference on the instabil i ty stri p of j3 Cep and SPB stars 
was studied bv iPamvatnvkhl il999T) . In that study, however, 
OP and OPAL opacity tables had also small differences in 
heavy elements abundances: Ni 5% more abundant in OP 
than in OPAL, and Ne and Fe ~ 2 — 3% more abundant in 
OP. 

Ni and Fe are the main contributors to the Z-bump, and 
Ni contributes to the Z-bump in kr a t higher temperatures 
than Fe (see e.g. lJefferv fe Said l2006l . In order to separate 
metal mixture effects from opacity calculation ones, we ob- 
tained OP opacity tables from the OPserver 1 for the same 
GN93 and AGS05+Ne metal mixture as OPAL. The abun- 
dances of the four elements that are not included in OP (P, 
CI, K, Ti) are redistributed among the other 15 elements. 
Doing so, the abundance differences between OP and OPAL 
opacity table computation are at maximum 0.06%. 

In Fig we show the difference in opacity for the inter- 
nal structure of two stellar models with the typical masses 
of an SPB (4 M ) and a (3 Cep (10 M ). We see that for 
the lower-mass model the opacity differences in the Z-bump 
region are of the order of 20%, and even larger for the 10 Mq 
model. 



4 STELLAR MODELS 

Stellar models have been computed with the code CLES 
(Code Liegeois d'Evolution Stellaire. IScuflaird I2005L The 
main physical inputs are : OP AL2001 equation o f state 
jRoeers fe Navfonovl EooS) and ICauehlan fe Fowlerl dl988l) 
nuclear reaction rates with lFonriicol£^^lTT200^r for the 
14 N(p,7) 15 cross-section. Convective transport is treated 
by using the classical Mixing Length Theory of convection 
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( Boh m-Vitensell958f) . and a convective overshooting param- 
eter of 0.2 pressure scale height was assumed in all the mod- 
els. For the chemical composition we have considered: GN93 
and AGS05+A(Ne)=8.11. We have computed models with: 
i) OPAL opacity tables with GN93 and it) AGS05+Ne chem- 
ical composition, then models with Hi) OP opacity tables as- 
suming GN93 and iv) AGS05+Ne mixtures. All the opacity 
tables are completed at logT < 4.1 with th e corresponding 
GN93 and AGS05 low temperature tables bv lFerguson et all 
(2005). The masses considered span from 2.5 to 12 Mq, and 
the chemical composition considered are: X = 0.70 for the 
hydrogen mass fraction and three different metal mass frac- 
tions, Z = 0.02, 0.01 and 0.005. For all the models the evo- 
lution was followed from the Pre-Main Sequence. 



5 STABILITY ANALYSIS 

Since oscillations in f3 Cep and SPB stars are driven by the 
K-mechanism acting in the Z opacity bump, the location of 
the instability strip in the HR diagram and the frequency of 
the excited modes are determined by the properties of the 
metal opacity bump. In models with the new OP opacity 
tables, the driving region, that is, the region where kt = 
(9 log nn/d log T) increases outwards, is found deeper in 
the star with respect to the models computed with OPAL. 
As a consequence, we expect hotter B-type pulsators with 
OP models compared to OPAL ones. 

We perform stability analyses of main-sequence mod- 
els from our grid using the non-adiabatic code MAD 
JPupret et al J 120031. Fig. H shows different effects of OP 
and OPAL on the p-modes excitation (upper panels) and 
on the high-order g-modes excitation (lower panels) for main 
sequence models of a 10 Mq star with initial chemical com- 
position given by X = 0.70, Z = 0.02 and GN93 metal 
mixture. For p-modes (or low-order g-modes), we see that 
in the common instability region, OPAL and OP predict the 
same excited modes, and their frequencies are not affected 
by opacity differences. There is, however, a significant effect 
on the hottest unstable models. While OPAL models with 
excited radial modes are expected only for X c < 0.5, ZAMS 
OP models have I — modes excited. 

The differences between OP and OPAL are even more 
pronounced for high-order g-modes. There are many more 
excited modes for OP models than for the OPAL ones. For 
instance, for a 10 Mq star, the OP model with X c = 0.2 
presents excited modes of i — 2 going from gie to (721, and 
the TAMS model has I = 1, <73i-<;43 and I = 2, <?28 _ <?53 
excited. In OPAL models, instead, no I = 1 high-order g- 
mode is excited; the first excited I = 2 g-modes appear at 
A c = 0.1, and in the TAMS model the order of the excited 
g-modes goes from n = 34 to n — 48. The T e g domain 
for which we find excited SPB-modes is therefore ~3000 K 
larger than for OPAL models. As a consequence, the number 
of expected hybrid /3 Cep-SPB objects is also larger for OP 
models (see right panels in Fig. [5] and Fig UK and c). These 
results are similar to those obtained bv IPamvatnvkhl (I1999T) 
with the previous OP opacity tables. In fact the higher Fe 
and Ni abundances in the original OP mixture, are balanced 
now b y a 18% higher opacity at the Z-bump. IPamvatnvkhl 
(1999) noted that the shift of the Z-bump towards higher 
temperatures leads to an increase of the /3 Cephei and SPB 
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Figure 2. Frequencies of pulsation modes as a function of logT e ff for main-sequence models of a 10 Mq Z=0.02 star. Unstable modes 
are described by large black dots and gray squares when models are computed, respectively, with OPAL-GN93 and OP-GN93 opacity 
tables. Lower panels represent high-order g-modcs. 
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Figure 3. As in Fig.[3]but comparing models computed with OPAL-GN93 and OP-AGS05+Ne opacity tables. 



instability domains towards higher T e g and luminosity with 
respect to OPAL models. These results, however, have not 
been taken into consideration to explain the high order g- 
modes detected in the /3 Cep stars 12 Lac and v Eri. 

The instability strips for OP and OPAL models, with 
GN93 mixture and metal contents Z=0.02 and Z=0.01 are 
shown in Fig. 2] (panels a and b). We see that the impact 
of OP is more important for lower metallicity. In fact, the 
expected OP (3 Cep pulsators in the LMC (Z=0.01) would 
be more than twice larger than estimated from OPAL stellar 
models. 

The Fe-mass fraction enhancement in the AGS05+Ne 
mixture, compared with GN93, has the main effect of ex- 



tending towards higher overtones the range of excited fre- 
quencies. Fig.|^]shows the t = Q,l and 2 excited p-modes for 
10 Mq main sequence models calculated with OPAL opac- 
ity tables for the two different metal mixtures. AGS05+Ne 
models also provides slightly wider instability bands, and 
this effect increases as metallicity decreases (see Fig. |1J3 and 
d), thus, the expected (3 Cep pulsators for LMC predicted 
with AGS05+Ne would be more than three times larger than 
with GN93. 

Finally, if we consider the combined effect of OP opac- 
ities and AGS05+Ne metal mixture, we see that the large 
difference in kr described in Fig. Q is reflected in wider in- 
stability strips, in the excitation of higher overtones in /3 Cep 
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Figure 4. Instability strips represented in a log T c ff-log P dia- 
gram. In each panel, the two regions of unstable modes represent 
j8 Cep- and SPB-type pulsations. 



models and in a significantly larger number of /3 Cep pul- 
sators for Z=0.01 (see Fig.0^ and f). 

Computations for the lower metallicity corresponding to 
SMC (Z=0.005), show that none of the different OP/OPAL 
and GN93/AGS05+Ne evolutionary tracks for masses up to 
12M0 predicts (3 Cep pulsators, whereas we find SPB-type 
modes excited when considering OP with AGS05+Ne. 



6 CONCLUSIONS 

We have shown that the current uncertainties on opacity cal- 
culations and on the solar metal mixture adopted in stan- 
dard stellar models, have a considerable effect on the ex- 
citation of pulsations in ft Cep and SPB stars. Compared 
to models computed with OPAL opacities and GN93 metal 
mixture, we find that with OP opacities high-order g-modes 
are predicted to be excited in hotter stars and that higher 
overtones are excited if the AGS05+Ne metal mixture is 
considered. For low metallicities, such as Z=0.01, a larger 
number of models with excited modes is also found. These 
findings could help solving, at least partly, the discrepancies 
between theoretical predictions and observations of (3 Cep 
in low-metallicity environments and of a large domain of ex- 
cited modes detected in some j3 Cep and hybrid SPB-/3 Cep 
pulsators. The detailed modelling of single stars, such as v 



Eri is, however, beyond the scope of this letter and will be 
addressed in a forthcoming paper. 
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